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SIZES AND SURFACE BRIGHTNESS PROFILES OF QUIESCENT GALAXIES AT Z - 2. 
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ABSTRACT 

We use deep Hubble Space Telescope Wide Field Camera 3 near-infrared imaging obtained of the GOODS- 
South field as part of the CANDELS survey to investigate a stellar mass-limited sample of quiescent galaxies at 
1 .5 < z < 2.5. We measure surface brightness profiles for these galaxies using a method that properly measures 
low surface brightness flux at large radii. We find that quiescent galaxies at z ~ 2 very closely follow Sersic 
profiles, with n mec u an = 3.7, and have no excess flux at large radii. Their effective radii are a factor ~ 4 smaller 
than those of low-redshift quiescent galaxies of similar mass. However, there is significant spread in sizes 
(ciog 10 r e = 0.24), with the largest z ~ 2 galaxies lying close to the z = mass-size relation. We compare the 
stellar mass surface density profiles with those of massive elliptical galaxies in the Virgo cluster and confirm 
that most of the mass-growth which occurs between z ~ 2 and z = must be due to accretion of material onto 
the outer regions of the galaxies. Additionally, we investigate the evolution in the size distribution of massive 
quiescent galaxies. We find that the minimum size growth required for z ~ 2 quiescent galaxies to fall within 
the z = size distribution is a factor ~ 2 smaller than the total median size growth between z ~ 2 and z = 0. 



Subject headings: cosmology: observations 
redshift 



1. INTRODUCTION 



galaxies: evolution — galaxies: formation — galaxies: high- 



Quiescent galaxies make up a consider able fraction of the 
massive galaxy population at z = 2 (e.g., iFranx et ail 120031; 
iDaddi et al.ll2005t iKriek et al.ll2006l) . Their structural evolu- 
tion has been the subject of considerable discussion, focus- 
ing in particular on their extremely compact natur e compared 



to low redshift galaxies of similar mass (e.g., | Daddi et al 



2005 



2008 



van Dokk um et aLM 2008: 



Hop kins et aflT2009t ISaracco et al.ll2009t 
2009b: iCassata et al.ll2010HMancini et alj|2010t iCassata et ail 
201 ll) . The early formation and subsequent evolution of these 
massive, compact objects presents a considerable challenge 
to current models of galaxy form ation and evolution (e.g., 
iWuvts etal.ll2010t I0seretalj|2011l) . It is unclear what the 
structure of the progenitors of these galaxies is, and the lack 
of extremely compact massive galaxies at low redshift im- 
plies considerable size evolution betw een z = 2 and z = 
dTruiillo et aLll2009b iTavlor et al J 120101) . However, efforts to 
accurately quantify this evolution are hindered by uncertain- 
ties. The apparent compactness of z ~ 2 quiescent galaxies 
may simply be an observational effect: photometric masses 
may be systematically overestimated due to modeling uncer- 
tainties, and siz es may be underestimated due to a l ack of 
imaging depth (Ho pkins et al.ll2009l:lMuzzin et al.ll2009l) . 

Due to the difficulty of obtaining high-quality spectra of 
quiescent galaxies at z > 1.5, dynamical m asses have only 
been measured for a few such galaxies ( Cappella ri et al.1 
20091 ICenarro & Truiilld 120091 Ivan Dokkum et a!T T2009b; 
Onodera et al.H2.Q10t Ivan de Sande et al.112.01 lb . Instead, pho- 
tometric stellar masses are used, which are subject to con- 
siderable uncertainties due to e.g., the quality of the stel- 
lar libraries used in modeling the spectral energy distribution 
(SED), or incorrect assumptions about the shape of the initial 
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mass function (IMF). These uncert ainties can result in sys- 
tematic errors of up to a factor ~ 6 dConrov et al.ll2009h . At 
low redshift there is good agreement between stellar masses 
determined b y photometric SEP fitting methods and dynam- 
ical masses (ITavlor et al.l 120 ldl) . Whether this is also the 
case at high redshift is unclear (e.g..lvan de Sande et al.ll201 ll 
Bezanso n etal .1120 1 It iMartinez-Manso et al.ll201 ll) . 

The second large source of uncertainty lies in the size de- 
termination of these galaxies. The compact objects observed 
at z ~ 2 may be surrounded by faint extended envelopes of 
material, which could be undetected by all but the deepest 
data. Stacking studies have been used to obtain constraints 
on the av erage surface brig htness profile of compact galax- 
ies (e.g.. Ivan der Wei et aTJ 12001 Ivan Dokkum et al.1 l2008t 
Cassa ta et al.ll2010b . However, detailed analysis of individual 
galaxies is more difficult, primarily due to the limited number 
of compact galaxi es for which ul t radeep near-infrared (NIR) 
data are available. ISzomoru et al.1 (120101) carried out an analy- 
sis on a z = 1 .91 compact quiescent galaxy in the Hubble Ultra 
Deep Field (HUDF) and confirmed its small size; 

In t his Paper we expand the analysis of ISzomoru et al.l 
(2010) using a stellar mass-limited sample of 21 quiescent 
galaxies. We make use of deep Hubble Space Telescope 
Wide Field Camera 3 (HST WFC3) data from the CANDELS 
GOODS-South observations to investigate the surface bright- 
ness profiles of quiescent galaxies at z ~ 2. These observa- 
tions are not as deep as the HUDF data, but cover a much 
larger area, allowing us to study a statistically more meaning- 
ful sample. We measure the surface brightness profile of each 
individual galaxy and investigate deviations from Sersic pro- 
files. Additionally, we compare the size distribution and pro- 
file shapes of z ~ 2 galaxies to those of low redshift quiescent 
galaxies. Throughout the Paper, we assume a ACDM cosmol- 
ogy with Q m = 0.3, fl A = 0.7 and H = 70 km s" 1 Mp c" 1 . All 
stella r masses are derived assuming a Kroupa IMF (iRroupal 
2001). All effective radii are circularized and magnitudes are 
in the AB system. 
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FIG. 1. — Left panel: rest-frame U — V and V — J colors of galaxies in the CANDELS GOODS-South deep field at 1.5 <£ < 2.5. Right panel: specific star 
formation rates as a function of redshift. Arrows indicate upper limits. The dashed line indicates where the specific star formation rate is equal to 0.3/?//. 
Quiescent galaxies selected using the UVJ color criterion are shown as filled green circles. Galaxies which are selected as quiescent based on their SSFRs are 
shown as open blue circles. There is good agreement between the two selection criteria. Both the J/VJ-selected galaxies and the SSFR-selected galaxies are 
included in our quiescent galaxy sample. 



2. DATA AND SAMPLE SELECTION 

We use NIR data t aken with HST WFC3 as part of 
the C ANDELS survey dGrogin etal.H20lU IKoekemoer et all 
|20T1 . This survey will target approximately 700 square ar- 
cminutes to 2 orbit depth in F105, /125 and H^q (COSMOS, 
EGS and UDS fields), as well as ~ 120 square arcminutes 
to 12 orbit depth (GOODS-South and GOODS-North fields). 
These NIR observations are complemented with parallel HST 
ACS exposures in VI506 an d /si4- We use the d e epest publicly 
available data, reduced by Koeke moer et al .1 (1201 ll) . which 
consist of /gi4, /125 and Hk,o observations to 4-orbit depth of 
a ~ 60 square arcminute section of the GOODS-South field. 
The full width at half-maximum of the point-spread func- 
tion (PSF) is w 0.18 arcsec for the WFC3 observations and 
« 0.11 arcsec for the ACS observations. The images have 
been drizzled to a pixel size of 0.06 arcsec for the WFC3 
observations and 0.03 a rcsec for the ACS observations (see 
IKoekemoer et all d201 ll) for details). 

Galaxies are selected in the GO ODS-South field us ing the 
^-selected FIREWORKS catalog dWuvts et al.H2008l) . This 
catalog combines observations of the Chandra Deep Field 
South ranging from ground-based f/-band data to Spitzer 24 
data, and includes spectroscopic redshifts where avail- 
able, as well as photo metric redshifts derived using EAZY 
dBrammer et al.1 120081) . These photometric redshifts have a 
median Az/(1 +z) = -0.001 with a normalized medi an abso- 
lute deviation of ct NM ad = 0.032 dWuvts et al.ll2008l) . Stellar 
masses were estimated from SED fits to the full photometric 
data set (N. M. Forster Schreiber et al. 201 1, in preparation), 
assuming a Kroupa IMF a nd the stellar population models of 
iBruzual & Criarlo1(l200"3l) . 

We select all galaxies with 1 .5 < z < 2.5 and stellar masses 
above 5 x 10 10 M(^, which is the completeness limit in this 
redshift range (Wuv ts et al.1 12009 ). In order to ensure that 
we include all quiescent galaxies we explor e both a color- 
color selection (the UVJ selection described in Willia ms et al.1 



(2009)) and a selection based on specific star formation rate 
(SSFR). In the left panel of Figure [TJ we show the rest-frame 
U — V and V-J colors of all z ~ 2 galaxies in the field. 
The das hed lines indicate the q uiescent galaxy selection lim- 
its from Willia ms et al.1 (120091) . Galaxies which fall within 
the dashed lines (green dots) have SEDs that a re consistent 
with red, quiescent galaxies. Patel et al] d201 ll) have shown 
that this selection method is very effective at separating dust- 
reddened starforming galaxies from truly quiescent galaxies. 
As an alternative to the UVJ selection we also select galaxies 
based on their SSFR. In the right-hand panel of Figure [T] we 
show the SSFRs of galaxies as a function of redshift. The SS- 
F Rs are estimated from the UV and 24 /im fluxes, as discussed 
in IWuvts et al.1 (I2009I) . The dashed line shows our selection 
limit, below which the SSFR is lower than 0.3/?//, where tu 
is the Hubble time. There is generally very good agreement 
between the two selection criteria, although several galaxies 
that seem to be quiescent based on their SSFRs are not se- 
lected by the UVJ method, and vice versa. We find no sig- 
nificant difference in the distribution of structural parameters 
of galaxies selected by either method; the median values are 
equal to within 6 percent, for the effective radii, Sersic indices 
and axis ratios. This is expected, given the large overlap be- 
tween the two samples. Since we wish to be as complete as 
possible we combine the two selection methods and include 
all galaxies selected by either method. This results in a sample 
of 21 quiescent galaxies, whose properties are summarized in 
Table [1] 

To illustrate the effects of our selection on galaxy morphol- 
ogy we show color images of all galaxies with 1.5 < z < 
2.5 and M s tellar > 5 x 10 10 M Q in the stellar mass-size plane 
in Figure [2] The color images are constructed from PSF- 
matched rest-frame t/336, B438 and #475 images, obtained by 
interpolating between the observed /gi4, /125 and H\^q im- 
ages. Although we do not select based on morphology, the 
galaxies in our quiescent sample (indicated with red crosses) 
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FIG. 2. — Stellar masses and sizes of galaxies at 1.5 < z < 2.5 with M ste u ar > 5 X 10 Mq. Color images are composed of rest-frame f/336, £438 and £475 
images, obtained from observed /814, /125 and H^q images. Galaxies which are included in our quiescent sample are indicated with red crosses. Although we do 
not select based on morphology, almost all galaxies in our quiescent sample are compact, bulge-dominated, and have red colors. 



are generally very compact, bulge-dominated systems with 
relatively red colors. Interestingly, all starforming systems 
at z ~ 2 appear t o have a well-defined red core, as was 
also pointed out by [ Szomoru et al. (2011) ( but also see, e.g ., 
iForster Schreiber et alj|201 Tat Forster Schre iber et alj|20lll) . 

3. MEASURING SURFACE BRIGHTNESS PROFILES 

Obtaining surface brightness profiles of high-redshift galax- 
ies is difficult, in large part due to the small size of these galax- 
ies compared to the PSF Direct deconvolution of the observed 
images is subject to large uncertainties. A common approach 
is therefore to fit two-dimensi onal models, convolved with a 
PSF, to the observed images. Isirsicl ( fl968l) profiles are com- 
monly used, since these have been shown to closely match 
the surface brightness profiles of nearby early-type galax- 



ies (e . g.,|C aon et HS1993L [Graham et al. 2003; Truiillo etalJ 
[2004tlFerrarese et alJl2006tlC6te et alJl2007tlKormendv et al. 
2009). However, there is no reason that high-redshift galaxies 
should exactly follow Sersic profiles. 

An obvious way to account for deviations from a Sersic 
profile is by using double-component fits, in which the devi- 
ations are approximated by a second Sersic profile. Although 
this provides a closer approximation to the true surface bright- 
ness profile than a one-component fit, it still depends on as- 
sumptions regarding the shape of the profile. We therefore 
use a technique which is more robust to deviations from the 
assumed model and accurately recove rs the true intrinsic pro - 
file. This technique was first used in Szo moru et alj ( 120101) : 
we summarize it here. First, we use the GALFIT package 
dPenp et alj|2002l) to perform a conventional two-dimensional 



TABLE 1 
Galaxy properties. 



ID a 




z 




R.A 


Dec 


stellar 

(logM Q ) 


SSFR 
(logyr 1 ) 


U- 


" Vrest 


V-J I!Xt 


(AB) 


(kpc) 


n c 


b/a 


P.A. d 

(deg) 


1060 


2 


.345* 


53 


.069829 


-27.880467 


11.14 


-10.16 


1 


.75 


0.92 


22.21 ±0.05 


2.75 ± 1.60 


9.21 ±1.10 


0.70 ±0.01 


-45.2 ±1.6 


1088 


1 


.752* 


53 


.065570 


-27.878805 


10.75 


-10.12 


1 


.28 


0.94 


21.84±0.03 


0.83±0.11 


5.50 ±0.67 


0.87 ±0.02 


-62.5 ±13.1 


1289 


1 


.759* 


53 


.116186 


-27.871904 


11.00 


-10.51 


2 


.04 


1.50 


22.35 ±0.02 


1.20 ±0.22 


3.26 ±0.40 


0.58 ±0.01 


-8.5 ±0.8 


1831 


1 


.536 


53 


.076366 


-27.848700 


11.25 


-9.42 


1 


.47 


1.21 


20.71 ±0.02 


2.34 ±0.33 


3.68±0.16 


0.92 ±0.01 


-44.6±3.2 


1971 


1 


.608 


53 


.150661 


-27.843604 


10.84 


-9.42 


1 


.63 


1.32 


21.71 ±0.03 


2.69 ±0.79 


5.07±0.31 


0.87 ±0.01 


27.0±2.1 


2227 


1 


.612 


53 


.150165 


-27.834522 


10.98 


-10.36 


1 


.54 


1.03 


21.40 ±0.02 


1.92 ±0.26 


3.76±0.15 


0.84 ±0.01 


-6.5 ±2.1 


2514 


1 


.548* 


53 


.151413 


-27.825886 


10.79 


-10.24 


1 


.67 


1.38 


21.96 ±0.06 


1.28 ±0.29 


5.73 ±0.93 


0.86 ±0.03 


7.0 ±5.9 


2531 


1 


.598* 


53 


.171735 


-27.825672 


10.87 


-10.51 


1 


.90 


1.28 


21.93 ±0.02 


1.18 ±0.12 


4.08 ±0.30 


0.95 ±0.02 


7.6± 11.6 


2856 


1 


.759* 


53 


.216633 


-27.814310 


10.83 


-10.37 


1 


.76 


1.54 


22.90 ±0.01 


0.87 ±0.03 


1.20 ±0.08 


0.63 ±0.02 


-16.5 ±1.0 


2993 


2 


.470 


53 


.163233 


-27.808962 


10.71 


-9.92 


2 


.35 


1.00 


23.38 ±0.02 


0.46 ±0.04 


1.01 ±0.09 


0.32 ±0.04 


-63.1 ±1.0 


3046 


2 


.125* 


53 


.116519 


-27.806731 


10.80 


-10.18 


1 


.51 


0.97 


22.22 ±0.02 


0.54 ±0.02 


3.59 ±0.34 


0.70 ±0.02 


-45.8±3.0 


3119 


2 


.349 


53 


.123107 


-27.803355 


10.94 


-9.85 


1 


.43 


0.82 


21.97±0.03 


0.47 ±0.06 


5.09 ±0.60 


0.49 ± 0.04 


79.2 ±1.5 


3242 


1 


.910 


53 


.158831 


-27.797119 


10.77 


-9.95 


1 


.45 


1.07 


22. 10 ±0.02 


0.39 ±0.03 


4. 17 ±0.45 


0.62 ±0.03 


57.1±3.0 


3548 


1 


.500* 


53 


.202356 


-27.785436 


10.76 


-10.36 


1 


.67 


1.37 


22.40 ±0.03 


0.67 ±0.05 


3.75 ±0.48 


0.65 ±0.04 


55.5 ±2.4 


3829 


1 


.924* 


53 


.069966 


-27.768143 


10.79 


-10.33 


1 


.90 


1.37 


22.85 ±0.05 


0.64±0.14 


4.24±1.15 


0.66 ±0.03 


44.7 ±3.2 


4850 


2 


.118* 


53 


.012891 


-27.705730 


11.17 


-10.36 


1 


.97 


1.67 


22.68 ±0.02 


0.84 ±0.09 


2.72 ±0.42 


0.20 ±0.02 


17.0±0.3 


5890 


1 


.756* 


53 


.174620 


-27.753362 


10.92 


-10.57 


1 


.76 


1.27 


22.09 ±0.01 


0.98 ±0.04 


1.89±0.12 


0.92 ±0.02 


-14.6±3.6 


6097 


1 


.903 


53 


.140997 


-27.766706 


11.21 


-10.48 


1 


.44 


1.22 


21.30 ±0.03 


1.86 ±0.43 


5.26 ±0.56 


0.79 ±0.02 


2.8 ±2.2 


6187 


1 


.610 


53 


.044923 


-27.774363 


11.71 


-10.58 


1 


.76 


1.87 


20.37 ±0.01 


7.08 ± 1.30 


2.77 ±0.05 


0.61 ±0.01 


-30.1 ±0.2 


6194 


1 


.605 


53 


.052217 


-27.774766 


11.19 


-10.68 


1 


.82 


1.47 


21.18±0.02 


1.97±0.14 


2.04 ±0.06 


0.57 ±0.01 


-56.3 ±0.4 


6246 


1 


.615 


53 


.043813 


-27.774666 


11.00 


-10.07 


1 


.76 


1.57 


21.71 ±0.03 


1.12±0.16 


10.10±2.38 


0.64 ±0.02 


-19.6±1.5 



"FIREWORKS ID (Wuyts et al. 2008) 

b Masses are corrected to account for the difference between the catalog magnitude and our measured magnitude. 

c Magnitudes, effective radii and Sersic indices are derived from the H[^q band residual-corrected profiles discussed in Section|4] 

d Position angles are measured counterclockwise with respect to North. 

*No spectroscopic redshifts are available for these galaxies; photometric redshifts are listed instead 
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Sersic profile fit to the observed image. For PSFs we use un- 
saturated stars brighter than = 22.86 that are not contami- 
nated by nearby sources. We verify the quality of our stellar 
PSFs by comparing their radial profiles to each other, and find 
that the profiles show small variations in half-light radius of 
order ~ 2%. We find no systematic dependence of these vari- 
ations with magnitude. In order to estimate the effects of PSF 
variations on our derived parameters we fit every galaxy us- 
ing each of the stars separately. We find that the derived total 
magnitudes, sizes and Sersic indices vary by about 0.1%, 3% 
and 7%, respectively. 

After fitting a Sersic model profile we measure the residual 
flux profile from the residual image, which is the difference 
between the observed image and the best-fit PSF-convolved 
model. This is done along concentric ellipses which follow 
the geometry of the best-fit Sersic model. The residual flux 
profile is then added to the best-fit Sersic profile, effectively 
providing a first-order correction to the profile at those loca- 
tions where the assumed model does not accurately describe 
the data. The effective radius is then calculated by integrating 
the residual-corrected profile out to a radius of approximately 
12 arcseconds (~ 100 kpc at z ~ 2). We note that the resid- 
ual flux profile is not deconvolved for PSF; however, we show 
below that this does not strongly affect the accuracy of this 
method. 

Errors in the sky background estimate are the dominant 
source of uncertainty when deriving surface brightness pro- 
files of faint galaxies to large radii. Using the wrong sky value 
can result in systematic effects. GALFIT provides an estimate 
of the sky background during fitting. To ensure that this es- 
timate is correct we inspect the residual flux profile of each 
galaxy at radii between 5 and 15 arcsec (approximately 40 to 
120 kpc at z = 2). Using this portion of the residual flux pro- 
file we derive a new sky value and adjust the intensity profile 
accordingly. We use the difference between the minimum and 
maximum values of the residual flux profile within this range 
of radii as an estimate of the uncertainty in the sky determina- 
tion. 

In ISzomoru et al.l d2010t) this procedure was tested us- 
ing simulated galaxies inserted into HST WFC3 data of the 
HUDF. Since the data used in this Paper are shallower we 
have performed new tests. We create images of simulated 
galaxies that consist of two components: one compact ellip- 
tical component and a larger, fainter component that ranges 
from disk-like to elliptical. The axis ratio and position angle 
of the second component are varied, as are its effective radius 
and total magnitude. The simulated galaxies are convolved 
with a PSF (obtained from the data) and are placed in empty 
areas of the observed Z/ieo band image. We then run the pro- 
cedure described above to extract surface brightness profiles 
and compare them to the input profiles. 

A selection of these simulated profiles is shown in Figure[3] 
The input profiles are shown as solid black lines. The dashed 
grey lines indicate the two subcomponents of each simulated 
galaxy. The directly measured profiles are shown in green. 
The best-fit Sersic models are shown in blue, and the residual- 
corrected profiles are shown in red. The residual-corrected 
profiles are plotted up to the radius where the uncertainty in 
the sky determination becomes significant. The effectiveness 
of the residual-correction method is clear: whereas a simple 
Sersic fit in many cases under- or overpredicts the flux at r > 5 
kpc, the residual-corrected profiles follow the input profiles 
extremely well up to the sky threshold (~ 10 kpc). The recov- 
ered flux within 10 kpc is on average 95% of the total input 



TABLE 2 
Surface brightness profiles. 



ID a 


'"arcsec 


' kpc 




logS 




(arcsec) 


(kpc) 


(AB mag arcsec -2 ) 


(logM kpc~ 2 ) 


1060 


0.0180 


0.147 


I8.413±0.0010 


10.843 ±0.0004 


1060 


0.0198 


0.162 


18.548 ±0.0011 


10.789 ±0.0005 


1060 


0.0216 


0.177 


18.673 ±0.0013 


10.739 ±0.0005 


1060 


0.0240 


0.196 


18.826 ±0.0015 


10.678 ± 0.0006 


1060 


0.0264 


0.216 


18.966 ±0.0019 


10.622 ±0.0007 


1060 


0.0288 


0.235 


19.095 ±0.0021 


10.570 ±0.0008 


1060 


0.0318 


0.260 


19.244 ± 0.0024 


10.510±0.0010 


1060 


0.0348 


0.285 


19.382 ±0.0027 


10.455 ±0.0011 


1060 


0.0384 


0.314 


19.534 ±0.0032 


10.395 ±0.0013 


1060 


0.0426 


0.348 


19.696 ±0.0037 


10.330±0.0015 



NOTE. — This Table is published in its entirety in the 
electronic edition of ApJ, and can also be downloaded from 
|http : / /www . strw . leidenuniv . nl/~szomoru/| A portion is 
shown here for guidance regarding its form and content. 
"FIREWORKS ID (Wuvts et al. 2008) 



flux, with a 1-(T spread of 2%. Recovered effective radii are 
less accurate, as this quantity depends quite strongly on the 
extrapolation of the surface brightness profile to radii beyond 
10 kpc. However, effective radii derived from the residual- 
corrected profiles are generally closer to the true effective 
radii than those derived from simple Sersic fits. 

4. MISSING FLUX IN COMPACT QUIESCENT Z - 2 
GALAXIES 

We now use the residual-correction method to derive the 
surface brightness profiles of the z ~ 2 quiescent galaxies. 
The results are shown in Figure @] The SEDs, shown in 
the top row, illustrate the low levels of UV and IR emis- 
sion of the quiescent galaxies in our sample. Rest-frame 
color images are shown in the second row. These images in- 
dicate that the galaxies in this sample generally have com- 
pact elliptical morphologies. Some galaxies have a nearby 
neighbor; in these cases we simultaneously fit both objects 
to account for possible contamination by flux from the com- 
panion object. In the third row, best-fit Sersic profiles are 
shown in blue and residual-corrected profiles in red. The 
residual-corrected profiles follow the Sersic profiles remark- 
ably well. Most galaxies deviate slightly at large radii. The 
difference between the best-fit Sersic profiles and the residual- 
corrected profiles are shown in the bottom row. The devia- 
tions are generally small within 2r e ; for some galaxies larger 
deviations occur at larger radii, but in these cases the un- 
certainty is very high due to the uncertain sky. Overall, the 
profiles are consistent with simple Sersic profiles. The pro- 
files are given in Table [2] and can also be downloaded from 
|http : / /www . strw . leidenuniv . nl/~szomoru/| 

In order to investigate whether the profiles of z ~ 2 quies- 
cent galaxies deviate systematically from Sersic profiles we 
plot the difference between the best-fit Sersic profile and the 
residual-corrected flux profile in Figure [5] for all galaxies. 
Black lines indicate the deviation profiles of individual galax- 
ies, and their mean is indicated by the red line. The light 
red area shows the l-a spread around the mean. The mean 
profile is consistent with zero at all radii; the surface bright- 
ness profiles of quiescent galaxies at z ~ 2 seem to be well 
described by Sersic profiles. On average the residual cor- 
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FIG. 3. — Effectiveness of the residual-correction for recovering surface brightness profiles. The method was tested on a large number of simulated galaxies, 
composed of two components: one compact bright component, and an extended fainter component. A small selection is shown here. The input profiles are shown 
in black, with the dashed grey lines indicating the two subcomponents. The PSF-convolved "observed" profiles are shown in green. Direct Sersic fits are shown 
in blue, and the residual-corrected profiles are overplotted in red. The shaded light red regions indicate the l-a errors due to uncertainty in the sky estimation. 
The size of the PSF half width at half maximum (HWHM) is indicated on the top axis of each panel. Input effective radii are indicated in black on the bottom 
axes. Effective radii derived from the direct Sersic fits and from the residual-corrected profiles are indicated in blue and red, respectively. The fraction of the 
input flux within 10 kpc recovered by the Sersic fits fseraic/^uiput i s given in each panel. The residual-corrected profiles clearly reproduce the input profiles more 
accurately than the simple Sersic fits, especially at large radii. 



rection increases or decreases the total flux of each galaxy 
in our sample by only a few percent, with an upper limit of 
7%. The mean contribution of the residual flux to the total 
flux for all galaxies in our sample is -0.7%. Thus, we do not 
find evidence that indicates that there is missing low surface 
brightness emission around compact quiescent z ~ 2 galaxies, 
and we therefore conclude that the small sizes found for these 
galaxies are correct. 

5. THE MASS GROWTH OF Z - 2 QUIESCENT 
GALAXIES 

In the previous Section we have shown that the surface 
brightness profiles of z ~ 2 quiescent galaxies closely fol- 
low Sersic profiles, and that their sizes are not systemati- 
cally underestimated due to a lack of sensitivity. We now 
compare their size distribution and surface brightness profiles 
to those of low-redshift galaxies. In Figure [6] we show the 
mass-size and magnitude-size relations for the z ~ 2 galax- 
ies a nd for low -redshift massive elliptical galaxies, taken 
from [Shen et al] ( 12003b . The z ~ 2 sample has been split 
into two redshift bins: 1.75 < z < 2.5 and 1.5 < z < 1.75 
(shown in blue and green, respectively). The low-redshift 



sample is shown in grey. Galaxies at z ~ 2 are significantly 
smaller than those at z = 0. We fit a power law of the form 
r e oc (l + z) ° and find a = -0 . 94 ±0 .16, which is compara- 
ble to e.g., Ivan der Wei et alJ (12008b and Ivan de Sande et ail 
(120111) . but slightly steepe r than iNewman et al l (120101) and 
significantly shallower than Buitrag o et al.1 d2008). 

However, the z ~ 2 galaxies span a large range in size; 
some are supercom pact, while other s are as large as z = 
galaxies. Following lShen et aD (120031) . we quantify this range 
using o"i g . , which is defined as the l-cr spread in log 10 r e 
around the median mass-size relation, which we fix to the 
z = slope. Note that we define th e scatter in logio basis, 
not the natural logarithm as used bv lShen et al.1 ( 12003b. It is 
equal to 0.24 ± 0.06 for our entire sample, while IShen et alJ 
(2003) find values around cx\ oe re = 0.16 for early-type galax- 
ies at z = 0.1 in the same mass range. The values for the two 
high-redshift subsamples are 0.21 ±0.11 at 1.5 < z < 1.75 
and 0. 19 ± 0.07 at 1 .75 < z < 2.5. These values are upper lim- 
its, since they include the errors on individual size measure- 
ments; however, if our error estimates are correct, their effect 
on the scatter is < 0.01 dex. The scatter we m easure is com- 
parable to that found in lNewman et alJ (1201 ll) . These authors 
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FIG . 4. — Broadband SEDs, color images and PSF-corrected surface brightness profiles of z ~ 2 quiescent galaxies. The SEDs, obtained with FAST I Kr iek et"al] 
2009), are based on photometry from the FIREWORKS catalog. The color images are composed of rest-frame Um, B438 and §475 images, obtained from the 
observed /814, /125 and //160 data. The red ellipses are constructed from the best-fitting effective radii, axis ratios, and position angles. The best-fit Sersic profiles, 
obtained using GALFIT, are indicated by blue dotted curves. Residual-corrected surface brightness profiles are shown in red. Effective radii and the PSF HWHM 
are indicated at the bottom and top axes, respectively. We are able to measure the true surface brightness profiles of these galaxies down to approximately 26 
mag arcsec~ 2 and out to r fa 10 kpc. In the bottom row we show the difference between the best-fit Sersic profile and the residual-corrected profile. Individual 
residual-corrected profiles show deviations from simple Sersic profiles, although these deviations are consistent with zero within the errors. 



find cr log|o ,. f w 0.25 for galaxies with 1O 1O,7 M < M steUar < 
10 1L7 M© at z ~ 2. We note that our sample contains sev- 
eral galaxies that are part of an overdensity at z = 1.6 ( e.g., 
iGilli et al.ll2003t lCastellano"etd]l2q07l: iKurk et al.H200l . In 
particular, the two largest galaxies in our sample are part of 
this overdensity. Excluding the z = 1 .6 galaxies from our anal- 
ysis does not significantly alter the spread in galaxy sizes in 
the 1.5 < z < 1.75 redshiftbin: ow n r . = 0-21 ±0. 14. 
The size measurements used in S hen et al.l {2003) have been 



shown to suffer from systematic errors due to background 
oversubtraction dGuo et al.l l2Q09f). As a result o f this, the 
mass-size relation measured by Ishen et alj (|2003l) is signifi- 
cantly shallower than that found by. e.g.. iGuo et alj (2009). 
We therefore repeat our determinatio n of the scatter a round 
the z ~ 2 mass-size relation using the lGuo et alj {2009) mea- 
surements. This results in a decrease in the scatter by only 
~ 0.03 dex, and does not affect our conclusions. 

We note that, even within the limited redshift range un- 
der consideration, differences in redshift play a role: the 





galaxies in the 1.75 < z < 2.5 subsample are clearly smaller 
than the 1.5 < z < 1.75 galaxies. This may explain some 
of the disagreement between studies of high-redshift quies- 
cen t galaxies. In particu lar, the large effective radii found 
by iMancini et alJ 42010) for some high-redshift quiescent 
galaxies could be due to the fact that they select galaxies 
with 1.4 < z < 1.75. In this context, part of the size evo- 
lution between z ~ 2 and z = could be due to the ap- 
pearance of yo ung, relatively large qui e scent galaxies af- 
ter z ~ 2 (e.g.. Ivan Dokkum et al.N2008t iFranx et alJ 12008b 
Saracco et alTl2009h Ivan der Wei et alJ l2009t ICassata et alJ 
2011b . Wenote that lSaracco et al.1 (120091) find evidence for 
a correlation of galaxy compactness with stellar age, such 
that the most compact high-redshift quiescent galaxies con- 



tain older stellar populations than quiescent galaxies that lie 
close to the z = mass-size relation. We investigate this corre- 
lation in Figure [7] using rest-frame U — V color as a proxy for 
galaxy age. We define galaxy compactness as the offset be- 
t ween the z ~ 2 gala xy sizes and the z = mass-size relation 
of lShen et al.1 (12001 : r e /r e . z=(l = r e /(2.88 x 10" 6 x M° 56 ). We 
find no evidence for a correlation between galaxy compact- 
ness and galaxy age in our data. 

In Figure [8] we compare the stellar mass surface density 
profiles of the z ~ 2 galaxies to those of low redshift galax- 
ies. Based on their masses and number densities, we expect 
z ~ 2 quiescent galaxie s to evolve into the most massive low- 
redshift galaxies (e.g., Ivan Dokkum et al1l2010t) . As a com- 
parison sample we therefore use surface brightness profiles 
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FIG. 5. — Deviations of galaxy profiles from Sersic profiles. The difference 
between the best-fit Sersic profile and the residual-corrected profile is plotted 
as a function of radius for all galaxies in our sample (black lines). The mean 
profile is shown in red, with the shaded light red region indicating the 1 — a 
spread in the distribution. Although individual galaxy profiles deviate from 
Sersic profiles, on average the difference is consistent with zero. 

of elliptical galaxies with equal or higher mass in the Virgo 
cluster from Korme ndv et al.l (12009b . These authors used a 
combination of space-based and ground-based observations 
to obtain surface brightness profiles with very high resolution 
and dynamic range, covering almost three orders of magni- 
tude in radius. The surface brightness profiles have been con- 
verted to stellar mass surface density profiles using the total 
stellar mass-to-light ratios. We have ignored radial color gra- 
dients, which are k n own t o exist at low and high re d shift (e.g., 
Ivan Dokkum et al.l (120101) : ISzomoru et ail d2011l) : iGuo et all 
(|2011|)). These profiles are shown in grey, with the profiles of 
the z ~ 2 galaxies overplotted in blue and green. 

What is most apparent in Figure [8] is that the central (r < 
1 — 3 kpc) surface densities of the z ~ 2 galaxies are very sim- 
ilar to those of the z = galaxies, while at larger physical radii 
(in kpc) the high-redshift galaxies have lower surface densi- 
ties than the low-redshift galaxies. Th e profiles are in close 
agreement with prev ious studies (e.g., Bez anson et al.ll2009l; 
iCarrasco et al1l20Toh . We compare the change in radial mass 
density pro files to the mass evolu tion of quiescent galaxies 
described in lBrammer et al.l (1201 11) . These authors show that 
galaxies with a number density of 10 -4 Mpc" 3 have grown in 
mass by a factor ~ 2 since z = 2. As mentioned above, the 
mass contained within 3 kpc changes very little from z ~ 2 
to z = 0; we find an increase on the order of 10%. However, 
the mass contained outside 3 kpc is approximately ten times 
higher for the z = galaxies than for the z ~ 2 galaxies, and is 
equal to 58% of their total mass. Thus, slightly more than half 
of the total mass of the z = ellipticals is located at r > 3 kpc, 
whereas the z ~ 2 galaxies contai n nearly no mass at th ese 
radii. This is consistent with the Brammer et al.l ( 1201 11) re- 
sult, and suggests that compact z ~ 2 quiescent galaxies may 
survive intact as the cores of present-day massive ellipticals, 
with the bulk of mass accretion since z ~ 2 occuring at large 
radii. This is consistent with a n inside-out scenario of g alaxy 
growth, as described in e.g., Ivan Dokkum et ail dMTol) . We 
note that this discussion ignores transformations of star form- 



ing galaxies to the quiescent population. 

Finally, we compare the comoving number densities and 
comoving cumulative number densities of our z ~ 2 sample 
to the number densities of z = galaxies in Figure [9] To ob- 
tain the z = number densities we c ombine the z = mass 
function for early-type g alaxies from iBell et al.l (120031) with 
the mass-size relation of IShen et aT] d2003l) : we use the re- 
lations appropriate for early-type galaxies and evaluate over 
the mass range 5 x 10 10 M Q < M stMar < 5 x 10 n M Q . Given 
our small field size we cannot determine number densities ac- 
cu rately. We therefore ad opt the number densities measured 
by Bramm er et al.l Q201 1). These authors used data covering 
a much larger field of view (approximately 25 times larger 
than the CANDELS GOODS-South field), and as such their 
results are less sensitive to cosmic variance. We scale our 
(cumulative) number density distri butions such that the t otal 
number density corresponds to the Brammer et al.l (1201 11) re- 
sults. We note that our measured number d ensities are ap- 
proxim ately a factor 2 smaller than those in iBrammer et al.l 
(120111). consistent with expec tations from field- to-field varia- 
tions (Somerville et al. 2004). We first consider the comov- 
ing number density distributions, plotted in the left panel of 
Figure [9] As expected, the median radius and the total num- 
ber density increase with time, as existing galaxies grow in 
size and new quiescent galaxies appear. r e me <n an = 0.84 ± 0.20 
kpc, 1.92 ±0.45 kpc, and 3.82 ±0.03 kpc at 1.75 < z < 2.5, 
1 .5 < z < 1 .75 and z = 0, respectively. 

We can place constraints on the minimum size growth of 
z ~ 2 galaxies by considering comoving cumulative number 
densities, shown in the right-hand panel of Figure [9] We as- 
sume that the population of z ~ 2 quiescent galaxies grows 
just enough to fall within the z = size distribution, but 
doesn't necessarily grow to the same median size as z = 0. 
This results in a shift of the z ~ 2 cumulative number den- 
sity distribution, indicated by the filled arrows in Figure [9] 
This shift is approximately a factor ~ 2 smaller than the size 
growth required for the z ~ 2 quiescent galaxy population to 
match the median size at z = (indicated by the open arrows). 
Thus, in this minimal-growth scenario, half of the observed 
size evolution between z ~ 2 and z = is due to the growth 
of existing galaxies, while the other half results from the ap- 
pearance of new, larger quiescent galaxies at int ermediate red- 
shifts. Th ese results are consis tent with e.g., iCassata et al.l 
d20Tl and lNewman et all ( l20TTb . 

6. SUMMARY AND CONCLUSIONS 

In this Paper we have demonstrated that the small measured 
sizes of z ~ 2 massive quiescent galaxies are not caused by a 
lack of sensitivity to low surface brightness flux. Using deep 
data and a method which is sensitive to excess emission at 
large radii, we have shown that the surface brightness profiles 
of these galaxies are well described by Sersic profiles. The 
median Sersic index is n me dian = 3.7, similar to low-redshift 
quiescent galaxies. 

The sizes of z ~ 2 quiescent galaxies span a large range; 
although the median effective radius is small (r eme di a „ = 1.1 
kpc), values up to ~ 7 kpc are observed. The scatter in log r e is 
0.24 at z ~ 2, aproximately 1 .5 times as large as at z = 0. This 
indicates that the "dead" population of galaxies is very diverse 
at z ~ 2. We note that the size evolution between z = 1 .5 and 
z = 2.5 is significant, which suggests that the cause of dis- 
crepancies in the results of different studies of the measured 
sizes of quiescent galaxies around z = 2 could be due to small 
differences in the redshift ranges considered. 
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FIG. 6. — Relations between size and stellar mass (left panel) and s ize and rest-frame r-band absolute magnitude (right panel). Grey lines indicate the low- 
redshift mass-size and magnitude-size relations from Shen et al. (2003), green and blue points indicate the z ~ 2 sample (divided into low and high redshift bins, 
respectively). The z. ~ 2 galaxies are, on average, almost an order of magnitude smaller than low-redshift galaxies of similar mass and luminosity. However, 
there is a significant range in sizes at both redshifts. The largest z ~ 2 galaxies lie very close to the z = mass-size relation. 
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FIG . 7 . — Offset from the z = mass-size relation as a function of rest-frame 
U — V color for galaxies at 1.5 < z < 1.75 (green points) and 1.75 < z < 2.5 
(blue points). The offset is calculated by dividing the effective radius of 
each galaxy by the median effective radius of z = quie scent galax i es wit h 
the same mass, using the z = mass-size relation from Shen et al. 1 2003). 
The galaxies are split into two redshift bins. Assuming that rest-frame 
U — V color is a good proxy for the mean stellar age of galaxies, we find no 
evidence for a correlation between galaxy compactness and galaxy age for 
z ~ 2 quiescent galaxies. 



FIG. 8. — Comparison of stellar mass surface density profiles of z ~ 2 
galaxies (blue and gTee n curves) to elliptical galaxies in the Virgo cluster 
(Kormendy et al. 2009; grey curves). The Virgo galaxies are selected to 
have masses equal to or higher than those of the high-redshift galaxies. 
Radial color gradients are ignored when calculating the mass density 
profiles. The star, top left, indicates the PSF HWHM at z = 2. The 
central densities of the z ~ 2 galaxies are very similar to those of the 
Z = galaxies. At larger radii, however, significant evolution must occur 
if the z ~ 2 galaxies are to evolve into massive low-redshift elliptical galaxies. 



Additionally, we have compared the stellar mass surface 
density profiles of z ~ 2 galaxies to those of massive early- 
type galaxies in the Virgo cluster. Although the densities 
within ~ 1 kpc are comparable, at larger radii the z ~ 2 galax- 
ies show a clear deficit of mass. This puts strong constraints 
on models of galaxy formation and evolution. Firstly, most of 
the size buildup of z ~ 2 quiescent galaxies must occur at large 
radii (> 1 kpc). Secondly, a significant contribution from ma- 



jor gas-rich mergers since z ~ 2 seems to be ruled out, as this 
would disturb the inner density profiles of these galaxies. Mi- 
nor, dry merging and slow accretion of matter seems to be the 
most viable method of evolving these galaxies into their z = 
descendants. 

Finally, we have investigated the evolution in the size dis- 
tribution of massive quiescent galaxies. We conclude that the 
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FIG. 9. — Comoving number density (left panel) and cumulative comoving number density (right panel) as a function of effective radius r e , for galaxies at 
1 .75 < z < 2.5, 1 .5 < z < 1 . 75, and z = (soli d blue, dashed green, and dot-dashed grey li nes, re spectively). The z = number densities are obtained by combining 
the stellar mass function of Bell et al. (2003) with the mass-size relation of Shen et al. (2003). The z 2 number densities have been scaled such that the total 
number density corresponds to the results of Brammer et al. 1 20 11). Both the median effective radius and the total number density of quiescent galaxies show a 
strong increase from z ~ 2 to z = 0. The solid arrows in the right-hand panel indicate the minimum size growth required for high-redshift galaxies to grow into 
the smallest galaxies at z = 0. The open arrows indicate the size growth required for high-redshift galaxies to grow to the same median size as galaxies at z - 0. 
The minimum size growth required for z. ~ 2 quiescent galaxies is approximately a factor 2 smaller than the median size growth between z ~ 2 and z. = 0. 



median size of massive quiescent galaxies changes by a fac- 
tor ~ 4 between z ~ 2 and z = 0, and is accompanied by an 
increase in number density of a factor ~ 7. However, it is 
important to note that the size growth of individual galaxies 
is likely to be significantly smaller. The minimum required 
size growth for the z ~ 2 quiescent galaxy population is ap- 
proximately a factor ~ 2 smaller than the median overall size 
growth. In this scenario the stronger overall size growth may 
be caused by the appearance of new, larger quiescent galaxies 
at intermediate redshifts. 

One of the main observational uncertainties pertaining to 
the size evolution of massive quiescent galaxies now appears 
to be resolved; robust sizes, measured at high resolution and 
using very deep rest-frame optical data, indicate that galaxies 
at z ~ 2 were significantly smaller than equally massive galax- 
ies at z = 0. However, the mechanisms driving this evolution 
and their precise effects on the structure of individual galax- 
ies, as well as on the characteristics of the population as a 
whole, are still not entirely understood. Most studies seem 



to point towards gas -poor galaxy merging as the dominant 
growth process (e.g ., iBezanson et al.ll2009FlNaab et al.ll2009b 
Hopki nset al.ll2010b ; however, it is unclear whether this can 
account for all the observed size growth. A complicating fac- 
tor in such studies is that tracing the same group of galaxies 
across cosmic time is very difficult, since their masses, sizes 
and stellar population properties are not constant; selecting 
the same population of galaxies at different epochs is there- 
fore not trivial. Studies at fixed (cumulative) number density 
may provide a solution to this problem, though only for rela- 
tively massive galaxies. 



This work is based on observations taken by the CAN- 
DELS Multi-Cycle Treasury Program with the NASA/ESA 
HST, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS5- 
26555. 

Facilities: HST (ACS, WFC3). 
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